Spin-orbit coupling induced Mott transition in Ca 2 _ x Sr a; Ru04 (0< x< 0.2) 



O 

(N 

C 

(N 



CO 

-4-* 



I 

C 

o 
o 



> 

-1- 

o 
in 

o 



Guo-Qiang Liu 

Max- Planck- Institut fur Festkorperforschung, D-70569 Stuttgart, Germany 
(Dated: June 27, 2011) 

We propose a new mechanism for the paramagnetic metal-insulator transition in the layered 
perovskite Ca2- a; Sr :E Ru04 (0<x<0.2). The LDA+!7 approach including spin-orbit coupling is used 
to calculate the electronic structures. In Ca2Ru04, we show that the spin-orbit effect is strongly 
enhanced by the Coulomb repulsion, which leads to an insulating phase. When Ca is substituted 
by Sr, the effective spin-orbit splitting is reduced due to the increasing bandwidth of the degenerate 
dxz and d yz orbitals. For x = 0.2, the compound is found to be metallic. We show that these results 
are in good agreement with the experimental phase diagram. 

PACS numbers: 71.30.+h, 71. 15. Mb, 71.27.H-a, 71.20.-b 



The layered perovskite Ca 2 -xSr x Ru04 (CSRO) has 
been intensely studied during recent years since this se- 
ries of compounds exhibits a variety of interesting physi- 
cal properties as a function of the Sr concentration li- 
Sr 2 Ru04 is a p-wave superconductor— £ with a K 2 NiF 4 - 
type structure. The substitution of Ca for Sr causes the 
RuOg octahedra to rotate, and start to tilt at x = 0.5£ 
Following with the structure distortion, CSRO under- 
goes a series of phase transition from a paramagnetic 
metal (0.5<x<2) to a magnetic metal (0.2<x<0.5), and 
finally to a Mott insulator (0<x<0.2)4 It is unusual that 
in the Mott insulating regime the metal-insulator tran- 
sition temperature (Tmi) is higher than the Neel tem- 
perature (Tjv) of the antiferromagnetic (AFM) phase, 
which shows that a paramagnetic (PM) insulating phase 
exists between these transition temperatures^ ' 10 ' 11 For 
pure Ca 2 Ru04, the PM insulating regime extends from 
T N = 110 K to T MI = 357 K^&ia This property makes 
Ca 2 - x Sr a; Ru04 (0<x<0.2) different from other AFM 
Mott insulators. 

Recently, Qi et al^ found that the substitution of the 
lighter Cr for the heavier Ru strongly depresses T M i in 
Ca 2 Rui_ 2/ Cr !; 04 (0<y<0.13), which implies a possible 
influence of the rclativistic spin-orbit (SO) coupling on 
the Mott transition as pointed out by the authors. It 
is well known that SO coupling plays an important role 
in 5d transition- metal oxides. For example, Kim et alr^ 
found that Sr 2 Ir04 is a J e ff = 1/2 Mott insulator, and 
they showed that the unusual insulating state can be ex- 
plained by the combined effect of the SO coupling and 
Coulomb interaction. In the 4d oxides, the importance 
of SO coupling is under debate. Mizokawa et a/.r* ob- 
served strong SO coupling in Ca 2 Ru04 from their pho- 
toemission experiment. Based this finding, they argued 
that the strong SO coupling in Ca 2 Ru04 would cause a 
complex electronic configuration. Theoretical studies re- 
vealed strong SO effects in Sr 2 Ru04 and Sr 2 Rh04r^^ 
which seemingly support the photoemission experiment. 
However, Fang et a/. 17 ' 18 reported an LDA+U study of 
Ca 2 Ru04. They found the AFM state has a rather sim- 
ple configuration xy ^ xz ^ yz ' without much influence 
of the SO coupling. These seemingly inconsistent view- 
points raise a question: what role does the SO coupling 



play in CSRO? 

In this paper we present electronic structure calcula- 
tions for Ca 2 _ x Sr a; Ru04 using the LDA+U method in- 
cluding the SO coupling. We show the combination of the 
SO coupling and Coulomb repulsion opens a band gap in 
PM Ca 2 Ru04. The appearance of the Mott insulating 
phase is strongly dependent on the tilting of the Ru06 
octahedra, which naturally explains the PM Mott tran- 
sition in the experimental phase diagram. On the other 
hand, we find SO has much less influence on the AFM 
order. We show that these phenomena can be explained 
by a simple formalism. 

All the calculations in this work were performed with 
the full-potential linear augmented plane wave (FLAPW) 
within the local-density approximation (LDA), as im- 
plemented in package WIEN2K^ Two experimental 
structure^ were considered in this work. For Ca 2 Ru04, 
we used the structure at 180 K, with the space group 
Pbca, lattice constant a=5.394, b=5.600, and c=11.765 
A£ For Cai 8Sro. 2 Ru04, we used the experimental struc- 
ture at 10 K, but the substitution of Sr for Ca is 
only taken into account via the structural changes. 
Ca1.8Sro.2R.uO4 also has the space group Pbca but with 
lattice constant a=5.330, b=5.319, and c=12.409 A£ For 
the AFM state, we considered the 'A-centered' mode. 3 
The LDA+U calculations were performed with U = 3.0 
eV, which is similar to the value used by Fang et o^ 17 ' 18 
We will show that this U value can reproduce the mea- 
sured band gap in Ca 2 Ru04. 

In Fig. [TJ we present our theoretical band structures 
for paramagnetic Ca 2 Ru04 using different approxima- 
tions. The LDA band structure is well know n 20 ' 21 : the 
bands crossing the Fermi level are from Ru i 2g orbitals, 
containing four d electrons. Our LDA band structure 
shown in Fig. la is consistent with the previous study^ 
The inclusion of the SO coupling (Fig. lb) only shows 
some slight changes on the band structure. This is not 
surprising since the SO coupling constant £ in Ca 2 Ru04 
is presumably similar to the one in Sr 2 Ru04, where it is 
only about 93 rneVJ^ The inclusion of Coulomb inter- 
action (Fig. lc) also shows little influence on the band 
structure since U does not break the orbital symmetry in 
the paramagnetic state. Surprisingly, the combined inter- 
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FIG. 1: Theoretical band structures for paramagnetic 
Ca2RuC>4 using different approximations, (a) LDA, (b) 
LDA+SO, (c) LDA+77, and (d) LDA+£/+SO. The LDA+U 
and LDA+C/+SO band structures are calculated with U= 3.0 
eV. 



action of the SO coupling and Coulomb repulsion gives 
a very different band structure compared to the LDA, 
LDA+SO or LDA+77 results. The LDA+C/+SO band 
structure shows an insulating phase with a gap about 0.2 
eV wide. The band gap obtained from the chosen U is in 
good agreement with the experimental data.— >22 

Similar combined effect of the SO coupling and U has 
been found in Si"2Rh04^ where it was termed Coulomb- 
enhanced spin-orbit splitting. Sr 2 Rh04 has a simi- 
lar crystal structure to Ca2_ x Sr x Ru04, and it can be 
regarded as a two-band (xz and yz) system since the 
xy band is below the Fermi level due to the RhOg 
rotation . 23 ! 24 The simpler problem of Sr2Rh04 can help 
us to understand the LDA+C/+SO band structure of 
Ca 2 Ru04. In Sr2Rh04, the SO coupling splits the de- 
generate xz and yz bands to the higher x±3/2 bands, and 
lower x±i/2 bands, where 

X3/2 = {xz + iyz) t, X_3/2 = (xz - iyz) I 
Xi/2 = {xz + iyz)l, X-i/2 = (xz - iyz) t • 



This splitting happens around the Fermi level, and there- 
fore the occupancies of the x±3/2 arL d x±i/2 states are 
changed: (n 1/2 + ^-1/2) - (%/2 + ^-3/2) = P > 0, where 



ni/2 



-1/2 



and 



n 3 /2 



-3/2- 



When the Coulomb 



interaction is taken into account, the SO splitting is en- 
hanced due to the different occupancies of the x±3/2 
and x±i/2 states. The interplay the SO coupling and 
Coulomb interaction can be represented by an effective 
SO constant^ 



teff = t+l(U-J)p, 



(1) 



where J is the Hund's coupling and p is determined self- 
consistently. 



> 1 

U 

c3 0.8 

O 06 

a 



^±1/2 " 
^±3/2 ■ 

xy 
2 2 
x -y 

o.2_ r 2 



-1.5 -1 -0.5 0.5 

Energy (eV) 



FIG. 2: Ru-d PDOS for paramagnetic Ca2Ru04 calculated 
by LDA+<7+SO. 



The problem of Ca2Ru04 is more complicated than 
Sr2Rh04 since the xy orbital is also involved. Fig. [2] 
presents the partial density of states (PDOS) for the Ru- 
d orbitals calculated by LDA+(7+SO. Here we present 
the PDOS for the x±3/2 and x±i/2 orbitals instead of 
xz and yz. The PDOS shows that the unoccupied t 2g 
bands (0.2-0.9 eV) are dominated by the x±3/2 states 
while the x±i/2 states are nearly fully occupied. The 
well separated x±3/2 and x±i/2 states indicate a large 
effective spin-orbit splitting in Ca2Ru04. Therefore, A 
simple explanation for the PM Mott transition is that 
the Coulomb-enhanced spin-orbit splitting opens a gap 
between the x±3/2 and x±i/2 bands, leading to an insu- 
lating phase with two holes residing on the x±3/2 orbitals. 
In this explanation, the xy state is assumed to be fully oc- 
cupied. However, in the experimental structure, the xy, 
xz and yz orbitals hybridize with each other due to the 
structural distortion. As may be seen, the weight of the 
xy state in the unoccupied t 2g bands is not small as shown 
in Fig. [2J The relative hole population shown in Fig. [5] 
is a;y:x± 1 :x±|=19 : 13:68, while this ratio is 21:39.5:39.5 
within LDA approximation. This shows that the inclu- 
sion of SO and U hardly changes the occupancy of the xy 
orbital. We may conclude that the band gap is mainly 
due to the splitting of x±3/2 and x±i/2 orbitals although 
the xy orbital is also involved in the Mott transition. 

Experimental research has found that the Mott tran- 
sition in CSRO is accompanied by an structural phase 
transition from the high temperature L-Pbca phase to 
the low temperature S-Pbca phase^ii where L (S) in- 
dicates a long (short) c-axis. The phase transition tem- 
perature Ts is a function of Sr concentration x, which 
decreases from 357 K at x = to K at x ~ 0.3^ii. 
For x > 0.2, CSRO is metallic and only has the L-Pbca 
phase. As indicated by Friedt et al.^, the structural tran- 
sition from the L-Pbca to S-Pbca phase is characterized 
by an increase in the tilting angle of RuOg octahedra. 
We will show that the tilting angle of RuOg plays an im- 
portant role in the Mott transition. To illuminate the 
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FIG. 3: Ru-t2 ff PDOS for paramagnetic state using different 
approximations. The left panels are for Cai.8Sro.2Ru04, and 
right panels for Ca2Ru04. 



relation between the Mott transition and the structural 
phase transition, we apply the LDA+U +SO calculation 
to Cai.8Sro.2Ru04. 

Fig. 3 presents our calculated PDOS for x = 0.2 (L- 
Pbca) and x = (S-Pbca). Fig. 3a and 3b show the LDA 
PDOS for x = 0.2 and x = 0. The xz/yz bandwidth is 
about 1.6 eV in Ca1.8Sro.2R.uO4, and it is reduced to 1.1 
eV in Ca2Ru04. The narrower xz/yz band in Ca2Ru04 
is due to its larger tilting angle. The tilting angle is about 
12° in Ca 2 Ru0 4 , and it is 6° in Cai. 8 Sr .2RuO 4 £ The 
tilting of the in-plane Ru-0 can significantly reduces the 



interaction between Ku-d xz / yz and 0-p z . Consequently, 
the xz/yz bandwidth decreases from x — 0.2 to x = 0, 
while the xy bandwidth is less influenced. With the nar- 
rower xz/yz band, Ca2Ru0 4 shows much higher xz/yz 
PDOS around the Fermi level than Cai.sSro.2Ru04. Fig. 
3c and 3d show the LDA+SO PDOS for x = 0.2 and 
x = 0. As may be seen, the occupancy difference be- 
tween the x±3/2 and x±i/2 states is larger in Ca2Ru04 
than in Cai.sSro.2Ru04. This is understandable if we 
consider the higher PDOS in Ca2Ru04. Eq. (1) shows 
the effective SO splitting is proportional to the occupancy 
difference p. Then the lager occupancy difference p in 
Ca2Ru04 will cause larger SO splitting when Coulomb 
interaction is taken into account. This is confirmed by 
the LDA+C/+SO PDOS shown in Fig. 3e and 3f. Using 
Eq. (1), we get Ce//=0.9 eV for Ca 2 Ru0 4 , and 0.6 eV 
for Cai.sSr .2RuO4. The larger SO splitting in Ca2Ru04 
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FIG. 4: Ru-t 2g PDOS for AFM Ca 2 Ru0 4 using different ap- 
proximations, (a) LDA+U, and (b) LDA+Zy+SO. 



leads an insulating phase, while Cai.gSro.2Ru04 remains 
metallic. Therefore, we have shown that the PM Mott 
transition in CSRO can be explained by the interplay of 
SO coupling, Coulomb interaction and structural distor- 
tion. 

CSRO is a single-layer system, where the t2 g bands 
are split into the singly degenerate xy band and doubly 
degenerate xz/yz bands. The bare SO coupling mainly 
influences the degenerate bands since Ru has a moderate 
SO constant. Therefore, the SO induced Mott transition 
is strongly orbital dependent. Fig. 3 indicates that the 
Mott transition in PM CSRO is driven by the narrowing 
of the xz and yz bands, while the xy orbital plays a 
less important role. This strong orbital-dependence could 
also be seen from Eq. (1). The bare SO parameter ( 
is constant for each orbital, but the Coulomb enhanced 
SO parameter is a function of orbital occupancies. 
This suggests that the strong orbital-dependence is an 
inherent feature of the SO induced Mott transition. 

Our calculation has shown that the PM insulating 
phase of CSRO originates in the strong effective SO split- 
ting. This picture supports the photoemission measure- 
ment by Mizokawa et al.— The suppression of T M i in 
Ca2Rui_ J/ Crj / 4 -i2 can also be understood within this 
picture. Cr has a much smaller atomic SO constant than 
Ru due to its smaller mass. And therefore the substitu- 
tion of Cr for Ru will reduce the SO splitting, leading to 
the observed decrease of T MI . 

As mentioned above, Fang et a/ i 17 ' 18 found that SO 
coupling has no much influence on the electronic con- 
figuration. They however pointed out that the pho- 
toemission measurement was done above the Neel tem- 
perature, while they applied the LDA+U method to 
the low temperature AFM state. To clarify if the 
SO coupling is less important in AFM state, we ap- 
ply the LDA+U and LDA+C/+SO calculation to AFM 
Ca2Ru04. Our LDA+U calculation gives a magnetic 
moment of mR U =l.25 ^b, which is consistent with Fang 
et al.'s calculation 17 ; while SO reduces the moment to 
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1.21 [is, showing a weak SO effect. The AFM PDOS 
for Ca 2 Ru04 are presented in Fig. 4. In contrast to 
the PM state, Fig. 4 shows that there is no Coulomb- 
enhanced SO splitting in the AFM state. The relative 
weak SO splitting in the AFM state can be explained by 
Eq. (1). The LDA+f/ calculation produces an insulat- 
ing phase for AFM Ca2Ru04 as shown in Fig. 4a. Since 
there is no density of states around the Fermi level, SO 
coupling can not change the orbital occupancies, which 
gives p = 0. Then we get £ e // = C> showing no enhanced 
SO splitting. It is noticeable that our calculations give 
an AFM ground state for Ca2Ru04, which is consistent 
with the experimental phase diagram. 5 

Comparing the PM and AFM insulating phases in 
CSRO, we may find that the two kinds of Mott transi- 
tion are similar. They both have an interaction to break 
the orbital symmetry. The interaction is SO coupling in 
the PM state and spin polarization in the AFM state. 
The breaking of orbital symmetry lifts the degenerate 
bands and changes the orbital occupancies. When the 
Coulomb interaction is taken into account, the orbital 
splitting, which is SO splitting in the PM state or ex- 
change splitting in the AFM state, is enhanced. If the 
enhanced splitting is large enough, it will lead to an in- 



sulating phase. The PM-AFM transition at T/v can be 
regarded as the competition of the Coulomb-enhanced 
SO splitting and the Coulomb-enhanced exchange split- 
ting. In the AFM state, the SO enhancement is quenched 
by the large exchange splitting, which causes the very dif- 
ferent electronic configuration from the PM state. 

In summary, we have applied LDA+L 7 +SO calcula- 
tions to CSRO. We find the Coulomb enhanced SO 
splitting produces an insulating phase in PM Ca2Ru04. 
This finding is consistent with the photoemission ex- 
periment, and also explains the recent experiment on 
Ca2Rui_ y Cr y 04. We show that the SO induced Mott 
transition in CSRO is driven by the change of the xz/yz 
bandwidth. For x = 0.2, the compound is found to 
be metallic. On the other hand, we find that SO cou- 
pling has much less influence on the AFM state, which 
is in agreement with the previous LDA+L/ study. The 
above picture shows that SO coupling plays a very sub- 
tle role in the correlated systems. The interplay of SO 
coupling, electron correlation and crystal structure dis- 
tortion would cause very rich physical phenomena. 

The author gratefully acknowledges Ove Jepsen for 
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